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Fast variability as a tracer of accretion regimes in black hole 
transients 
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ABSTRACT 

We present the rms-intensity diagram for black hole transients. Using observations taken with 
the Rossi X-ray timing explorer we study the relation between the root mean square (rms) am- 
plitude of the variability and the net count-rate during the 2002, 2004 and 2007 outbursts of 
the black hole X-ray binary GX 339-4. We find that the rms-fiux relation previously observed 
during the hard state in X-ray binaries does not hold for the other states, when different re- 
lations apply. These relations can be used as a good tracer of the different accretion regimes. 
We identify the hard, soft and intermediate states in the rms-intensity diagram. Transitions 
between the different states are seen to produce marked changes in the rms-fiux relation. We 
find that one single component is required to explain the ~ 40 per cent variability observed at 
low count rates, whereas no or very low variability is associated to the accretion-disc thermal 
component. 
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1 INTRODUCTION 

Black Hole X-ray transients (BHTs) are observed in different states 
during their outburst evolution. They correspond to particular spec- 
tral and timing properties seen in the observations. The physical 
interpretation of the states is still strongly discussed, but they are 
probably associated with different accretion regimes. A hard state 
[historically known as low/hard state (LHS)] is generally observed 
at the beginning and at the end of the outburst. It is roughly charac- 
terized by a power-law shaped energy spectrum with a photon index 
of ~ 1.6 (2-20 keV) and a high level of aperiodic variability [root 
mean square (rms) amplitude above ~ 30 per cent]. In the middle 
of the outburst, the energy spectrum is dominated by a soft thermal 
component and almost no variability is seen. During this soft state 
[historically known as high/soft state ( HSS)] a hard tail up to ~ 1 
MeV is also present dGrove et alj 1998). In contrast, a high-energy 
cut-o f f ( <; 200 keV ) is observed in the LHS (e.g. I Wilms et all 
2006; Motta et al. 2009). From a physical point of view, the ther- 
mal component present in the HSS is usually associated with emis- 
sion from an optically thick accretion disc, whereas the emission 
in the hard state is believed to arise from a 'corona' of hot elec- 
trons, where seed photons from an optically thin accretion disc are 
up-Comptonized. Synchrotron emission from a jet, which is known 
to dominate the radio and infrar ed spectrum during the LHS (e.g. 
lFendej2006l : lRussell et all2 006). co uld also have a signi ficant con- 
tribution to the hard X-ray emission jRussell et al1l2010h . 
In between these two 'canonical' states the situation is rather more 
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complex. Hard-to-soft and a soft-to-hard transitions are observed 
in relatively short time scales (hours/days) as compared to the ones 
seen for the canonical states (weeks/months). During these tran- 
sitions, both timing and spectral properties cha nge dramatica lly, 
leading to d ifferent states classifica tion s (see e.g. Belloni 20 id for 
a review). iHoman & Bellonll J2005t) and lBelloni et alJd2005l) (B05) 
identify two additional states, the hard-intermediate state (HIMS) 
and the soft-inter mediate state (SIMS) based on spectral and timing 
properties (see e.g Casella et al. 2004 for different types of QPOs). 
McCl intock & RemillarcT fcooe) propose a more quantitative but 
model dependent classification, which requires of a steep power- 
law state and an intermediat e state in ad dition to the two canonical 
states (for a comparison see lMotta et al.120091 ; hereafter M09). 
A fundamental tool for studying the evolution of BHT is the 
hardness-intensity diagram (HID), where the evolution of the spec- 
tral hard ness as a function of the acretion rate can be followed 
(see e.g. IHoman et al Jl200lh . The hard-soft-hard usual evolution 
and the main transitions become apparent in the HID. The latter 
are usually not associated with major changes in flux, with the fi- 
nal soft-to-hard transition occurring at count-rates ~ one order the 
magnitude lower. While the HID provides a general description of 
the BHT evolution, it is not enough for detailed studies (e.g. it is 
not able to establish with accuracy when the main transitions oc- 
cur). Complementary tools like power density spectra (PDS), the 
hardness-rms correlation (e.g. B05), and multiwavelength obser- 
vations are used for a more complete description of the different 
states and transitions. In this letter we present a new tool, the rms- 
intensity diagram, which allows one to map BHT states by only 
looking at the evolution of count-rate and rms (i.e. without spectral 
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Figure 1. Rms-intensity diagram for the 2002, 2004 and 2007 outbursts of GX 339-4 Dotted lines repres ent the 1, 5, 10, 20, 30 and 40 per cent fractional 
rms levels. The grey area marks the part of the diagram where solely type-B QPOs (Casella et al. 2004) are seen. The dashed line joins two consecutive 
observations separated by 33 days due to observability constraints. The cross corresponds to the first observation of the 2007 outburst. 



information). As a first step we have studied the case of the BHT 
GX 339-4, in which detailed observing campaigns have been per- 
formed during the multiple outburst observed with the Rossi X-ray 
timing explorer (RXTE). 

GX 339-4 is a low mass x-ray binary (LMXB) harbouring a > 
6 Ma accreting black hole ( Hynes et al. 2003 ; M unoz-Darias et all 
l2008h . Since its discovery ( IMarkert et all 19731) the system has un- 
dergone several outburst, becoming one of the most studied X-ray 
transient. The intense monitoring carried out by RXTE during the 
2002, 2004 and 2007 outburst has yielded detailed studies on the 
evolution of black hole states along the outburst (see e.g. B05). 
Here, we use this rich data set to study the evolution of the flux and 
variability along the outburst. 



2 OBSERVATIONS 

We have used all the RXTE public archival observations from the 
2002 (206), 2004 (295) and 2007 (239) outburst of GX 339-4. 
For the analysis we have considered only the data from the 
Proportional Counter Array (PCA). GOODXENON, EVENT and 
SINGLE-BIT data modes are used for the variability. Count-rates 
have been computed using PCA Standard 2 mode data correspond- 
ing to the PCU unit #2. The analysis has been performed using 



HEASOFT V. 6.7. and custom timing software running under IDL. 



3 THE RMS-INTENSITY DIAGRAM FOR THE 2007 
OUTBURST 

The rms-intensity diagram (RID) for the 2007 outburst of GX 339-4 
is presented in FigQ]as filled dots. For each observations net count 
rate corresponds to PCA channels 0-35 (2-15 keV). Only entire 
RXTE observations are initially considered, i.e. every dot corre- 
sponds to one observation. Power density spectra (PDS) for each 
o bservation have been also computed using the procedure outlined 
in lBelloni et alj l[2006j). We have used stretches 16 s long and the 
same energy selection as for the count rate. Total power was com- 
puted within the freq uency band 0.1-64 Hz and the fractional rms 
calculated following |Belloni_& Hasingei (1990). The absolute rms 
displayed in the diagram (FigQ} is obtained by multiplying frac- 
tional rms by net count-rate (PCU #2) for each observation. The 
black solid line joins observations contiguous in time starting from 
the observation marked with a cross (i.e. observation #1). We find 
that the source describe a continuous hysteresis pattern in the anti- 
clockwise direction. Four different regions can be distinguished in 
this diagram. 
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Figure 2. Same as Fig[T]but using fractional instead of absolute rms. Using 
this representation it is also possible to define a region where only type B 
QPOs are present. 

3.1 The hard line 

Starting from the first observation (cross in Fig[T) and following the 
solid line we see rms increasing with count rate. This linear trend is 
seen during the first 18 observa tions (i.e. ^ 40 day s). Following the 
state classification reported bv lMotta et alj J2009I) . all these obser- 
vations belong to the canonical LHS. In the context of this digram 
we will call this linear relation 'hard line' (HL). A continuous in- 
crease in hardness (M09) is observed during this HL. Weak, type-C 
QPOs appear in the PDS of the observations corresponding to the 
upper part of the HL without modifying the rms-flux relation. They 
become strong once the system shows fractional rms of ~ 30 per 
cent, and it is close to abandon the HL. In Fig. [2] we plot the RID 
but using fractional instead of absolute rms. In this representation 
the HL we see in Fig. Q] corresponds to the first 18 points starting 
from the cross, where fractional rms drops from 40 to 30 percent. 

3.2 State transitions 

From observation #19 to #23 (~ 3 days), the rms starts to decrease, 
whereas the flux is still increasing. Therefore, the system abandons 
the HL and follows a different linear relation. From observation #24 
to #29 (~ 3 days) rms a sudden increase in count rate is observed 
whereas rms is almost constant within ~4 counts s _1 . The spectral 
fits performed by M09 require of a thermal black-body component 
from observation #25 (see Fig. [I). Thus, we associate this increase 
in flux at constant rms to the appearance (2-15 keV band) of an 
optically thick accretion disc with a very low variability level. 

Until observation #32, type-C QPOs are observed in the cor- 
responding PDS. From observation #32 to #33 the rms decreases 
dramatically, crossing the 10 per cent fractional rms line. A type-B 
QPO is observed in the PDS. Although the change in rms is major, 
the observations are only separated by ~ 1 day. 

3.3 The soft branch 

From observation #34, the rms is below ~ 5 per cent and the count 
rate fades from ~ 1300 to ~ 100 cts s~ . A linear relation can be 
observed, but much more scattered than during the HL (see Fig.[]}. 
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Figure 3. Same as Fig[T]but using a log-log scale. The adjacent hard line 
(AHL) and the hard line at lower count rate becomes evident. The presence 
of a flux component with no variability is ruled out (see Sec. [6). Dotted lines 
represent the same fractional rms levels as in Fig[T] 



This new relation is followed by the system during the canonical 
soft state (M09). 

A sudden increase in rms is also observed during four observations. 
Three of these cross the ~ 10 per cent line and show type-C QPOs 
whereas a type-B is observed for the other. Observations with type- 
A QPOs are not well differentiated from the soft branch. However, 
some of them are placed close or slightly above the 5 per cent rms 
line, which roughly delimits the soft branch (see also Fig. [2}- 



3.4 Soft-to-hard transition. The adjacent hard line 

Fig. [3] shows the RID in log-log scale, where a monotonic rms in- 
crease is observed at constant flux (~ 100 cts s _1 ) from the bottom 
of the soft branch. Once the rms cross the ~ 20 per cent level a 
linear decrease in rms and flux towards the HL is observed. During 
this adjacent HL (see note in Fig. [3]) PDS are the typical of the 
LHS, although the spectral color is softer (i.e. hardness is lower) 
than the observed during the HL. 

4 THE 2002 AND 2004 OUTBURST 

In Figs. [Hill and [5] the RID corresponding to the 2002 and 2004 
outburst are shown as open triangles and open dots, respectively. 

• The 2002 RID is similar to the 2007 one (as in the case of the 
HID). The system leaves the HL roughly at the same position and 
an increase in flux at ~ constant (absolute) rms is also seen in the 
region where the ac cretion disc is observed by M09 during 2007. 
iHoman et iD d2005l) performed a multiwavelength campaign dur- 
ing the 2002 outburst. They find that the optical and near infrared 
properties of GX 339-4 change on MJD 52398, whic h correspond 
to the first observation outside of the HL (see Fig.fU- Belloni et al. 
(2005) also identify this observation as the first one belonging to 
the HIMS (i.e. non LHS). 

• The 2004 ourburst is fainter and therefore its absolute rms is 
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Figure 4. Soft (PCA 0-13; 2-6 keV) and hard (PCA 14-35; 6-15 keV) 
rms-intensity diagrams for the 2007 outburst of GX 339-4 plotted in log- 
log scale. The dotted lines represent the same fractional rms levels as in 
Fig[T] Since an accretion disc component is detected in the energy spec- 
trum (Motta et al. 2009; see arrows in the plot), the fractional rms level 
becomes higher in the hard band, the first observation showing a type-B 
QPO is indicated with a star. Crosses indicate the first observation. 

also significantly lower. In this case the variability does not de- 
crease when the system leaves the HL, but increases with flux fol- 
lowing a different linear relation. A clear feature related with the 
appearance of a thermal disc in the spectrum is not observed, sug- 
gesting a weaker disc component. Flux finally increases as rms de- 
creases until the systems reach the soft branch. After a long perma- 
nence in this branch rms increases up to ~ 20 per cent and goes 
back to the HL following the same adjacent HL observed in the 
other two outbursts. 

The HL is identical for the three outbursts, independent of the 
maximum flux reached. We also observe this behaviour in the BHT 
H1743-322 (Munoz-Darias et al. in prep). The soft branch is also 
very similar between the three outbursts, lying within the ~ 1 — 5 
per cent fractional rms region. We notice that observa tions just out- 
side t his region contains fast state transitions (see iNespoli et al.l 
l2003l ; ICasella et al]|2004l) . A detailed state classification in this re- 
gion of the diagram is beyond the scope of this work. 

It is also possible to define a region of the diagram (count rate 
> 400 counts s~ and ~ 7 — 10 per cent fractional rms) where 
only observations showing type-B QPOs (grey region in Fig. Q] 
and Fig. [2} are s een. These observations must be classified as 
SIMS according to lHoman & Belloni (2005). We also identify two 
observations (70110-01-47-00 and 70108-03-02-00) with fainter 
B QPOs at lower fractional rms than the region. A more detailed 
analysis of these data shows that the type-B QPO is only present in 
part of the observation, with rms being the expected for the QPO 
B region. 



5 HARD/SOFT RMS-INTENSITY DIAGRAMS 

With the aim of tracing the evolution of soft and hard variability 
we have computed the RID using soft (PCA 0-13; 2-6 keV) and 
hard (PCA 14-35; 6-15 keV) bands for the three outburst. These 



energy ranges can be selected independently of the observing mode 
in which data were taken, and allow us to cover the whole evolution 
of the three outbursts. In Fig.|4]the soft and hard RID are plotted for 
the 2007 outburst as open circles and open triangles, respectively. 
It is clear that both RIDs show the same HL, where fractional rms 
drops from 40 to 30 percent rms as the system becomes brighter. 
This is expected since is well-known that r ms is weakly energy 
depen dent (2-20 keV) during the LHS (see e.g. lMunoz-Darias et al.l 
l20ld for the cases of XTE J1752-223 and Cyg X-l). Differences 
arise when the system leaves the HL (observation #18; 92035-01- 
02-04): 

• During the next five observations the hard RID shows a de- 
crease in rms at constant count-rate. Then, once the disc compo- 
nent appears in the spectrum (M09; see arrow in Fig.0, a constant 
decrease in flux and rms following the 20 per cent line is seen until 
the type-B QPO appears (#33; star in Fig.[4}. For the same obser- 
vations an increase of flux at constant rms is observed in the soft 
RID. 

• Since the disc component is detected by M09 (see arrows in 
Fig. E]l , more fractional rms is seen systematically at higher ener- 
gies, even at the very bottom of the soft branch. 

On the other hand, the same adjacent HL is observed for both RID, 
as expected for LHS observations with flat rms spectra. 



5.1 The 2002 and 2004 outburst 

The behavior observed during the 2002 outburst (not shown here) 
is similar to the one described above. During 2004, the vertical rise 
after leaving the HL is much less evident in the soft RID, suggest- 
ing weaker accretion disc component. The decrease in rms at con- 
stant flux is neither observed in the hard RID, although the decay 
following the 20 per cent rms line is seen. 



6 DISCUSSION 

During the outburst evolution of BHT different accretions regimes 
are present. They yield the different states that we see when look- 
ing at the X-ray emission arising from these systems. In each one of 
these states, one or more physical mechanisms (e.g. thermal or non- 
thermal Comptonization in the corona, synchrotron emission from 
a jet, thermal emission from an accretion disc) play a role. The rel- 
ative contribution of these mechanisms can be studied using energy 
spectra, but also by analysing the particular variability they prompt 
in the light-curves. In this work we have shown that the evolution 
of the fast variability (0.1-64 Hz) can be used as a good tracer of 
the different accretion regimes in black hole binaries. This would 
be impossible in case the variability associated with the different 
spectral components were similar. Indeed, we find evidence for the 
presence of highly variable and an almost non variable components 
in the lightcurves of GX 339-4: 

• A sharp, linear relation between absolute rms and flux is 
observed at the beginning and at the end the three ourburst 
of GX 339-4 we have s tudied. This has been observed before 
by iGleissner et al. I J2004h in Cyg X-l during LHS observations. 
During these epochs, flux variations are probably tracing ac- 
cretion rate variations and therefore rms is probably correlated 
with accretion rate . Therms-flux relation was firstly observed by 
lUttlev & M cHardv < 200lh at much shorter time scales (~ 10 s) 
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than the ones presented here (~ 1 ks). As a test we have studied 
the short term variations of the rms-flux relation by dividing the 
lightcurve in 32 s segments for the first and the last HL observa- 
tions (2007 outburst). In both cases we also found a linear relation, 
with the fractional rms being consistent with the average of the en- 
tire observation. We note that a detailed study on these short-term 
variations is beyond the sc ope of this work. 

lUttlev & McHardvl d200lh found that the extrapolation of the rms- 
flux relation was consistent with no variability at lower count rates 
than the ones they studied. This supports the presence of two dif- 
ferent component in the lightcurves. GX 339-4 offers the opportu- 
nity to study this linear relatio n over a much broade r flux range. 
Consistently with the found bv lGleissner et alj d2004h by extrapo- 
lating the relation they obtain for Cyg X-l, we do not see remain- 
ing flux when the rms is close to zero (see Fig. [3}, thus our re- 
sults are consistent with the presence of one component solely re- 
sponsible for the variability at low count rates. Indeed, the system 
reaches its maximum variability level when this component dom- 
inates the spectrum . How ever, in the upper part of the HL, where 
lUttlev & McHardvl < l200ll) perform their study, the fractional rms 
decreases from 40 to 30 per cent. This could be associated with the 
presence of a second, less variable component. 

• We find strong evidence for the presence of an extra compo- 
nent associated with the appearance of thermal emission in the en- 
ergy spectrum of the source. This accretion disc component is not 
variable or varies very little. Its signature (i.e. flux increase at con- 
stant rms) is clearly seen in the 2002 and 2007 outburst once the 
systems leaves the HL. As expected, this becomes more evident if 
one only considers the 2-6 keV band (see Fig.[4j. Indeed, The frac- 
tional variability of the system is lower at lower energies from the 
first observation in which the thermal component is detected until 
it returns to the LHS at the end of the outburst. 

• The behaviour of the hard (6-14 keV) RID is complex. 
Once the system leaves the HL, a decrease in rms at constant 
flux is observed, suggesting that the properties of the physical 
mechanism that drives the hard emission are different. This could 
be explained by invoking a change in the properties of the corona 
or the base of the jet. For instance, the high energy cut-off is 
observed to increase during the hard-to-soft transition (M09), 
which could be interpreted in terms of a cooling of the corona (see 
M09 and reference therein for details). It is remarkable that very 
few observations show fractional rms < 5 per cent in the hard 
band, suggesting that an underlying variable component (probably 
associated to power-law emission) is present even during the soft 
state (see Fig.[4]l. 

The different rms-flux relations outlined above enable to per- 
form a state classification of the different observations using the 
rms-intensity diag ram. The following states can be identified using 
the description bv lHoman & BellonjHool) : 

• Hard line observations corresponds to the LHS, with rms in 
the range ~ 30-40 per cent. During these observations PDS are 
qualitatively similar between each other, and can be fitted by us- 
ing several broad Lorentzians (see e.g. Fig. 3 in iBelloni et aljiool 
where the three upper PDS correspond to observations at the begin- 
ning, the middle and the end of the HL during the 2002 outburst). 
We associated the sharp changes we see in the rms-flux relation at 
the end of the HL to the transition to and from the HIMS. This is 
supported by spectral and timing changes (e.g. M09, B05) and par- 
ticularly by the coincidence in time between the drop we see in r ms 
and that observed in optical and infrared bv lHoman etal] b005l) . 



• During the HIMS, the rms is observed to drop from ~ 30 to 
10 per cent. The presence of the accretion disc in the spectrum is 
detected when the fractional rms is ~ 20 per cent. We note that this 
point marks the transition between the hard and the steep power- 
law or intermediate states according to iMcClintock & Remillard 
d2006h (see M09). However, the major transition seems to occur 
before, when the system left the HL. 

• Soft-state observations are observed to follow a different rms- 
flux re lation, consistent with the results obtained bv lGleissner et al.1 
(2004) using relatively soft observations of Cyg X-l. The soft state 
in GX 339-4 is much softer than in Cyg X-l and a soft branch 
very different from the HL is clearly seen in the RID. It shows a 
variability level lower than ~ 5 per cent and much more scatter than 
the HL. This suggest fast changes in (at least) one of the spectral 
components of the lightcurve. 

• We find that between ~ 7 — 10 per cent fractional rms only 
type-B QPOs are observed (grey regions in Fig.[T]and|2](. They are 
SIMS observations according to lHoman & Bellonil J2005h and our 
diagram suggest that this kind of oscillations could follow a par- 
ticular rms-flux relation. We note that transitions between soft and 
SIMS are very fast, which result in hybrid observations where tim- 
ing properties are not constant. Thus, an accurate state classification 
close to the ~ 5 per cent rms line is in some cases not compatible 
with using entire 1-3 ks observations. 

• An adjacent hard line is observed at the end of the three out- 
bursts. Observations during this adjacent HL show PDS typical of 
the LHS but with a slighly softer spectrum than that of the HL. 
Why this other linear relation, not seen when the systems leaves 
the HL, is prese nt is unclear. It coul d suggest the presence of new 
components (see lRussell et alfeoioh or the total absence of any of 
the components (e.g. accretion disc) observed during the flux rise 
phase at the beginning of the outburst. 



7 CONCLUSIONS 

We have proved that the evolution of the fast variability can be used 
as a good tracer of the different accretion regimes in black hole 
binaries. In particular, we find that apart from the linear rms-flux 
found during LHS, different relations are followed by black hole 
binaries during the soft and intermediate states. 
In this work we have presented the rms-intensity diagram, showing 
that it is possible to associate the different regions of this diagram 
to the different states observed in black hole transients. Transitions 
can also be identified thanks to the marked changes they produce in 
the diagram. The application of this new tool to other systems and 
its interpretation will provide new insights in our understanding of 
the physical processes that take place during accretion episodes in 
X-ray binaries. 
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(FP7/2007-2013) under grant agreement number ITN 215212 
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